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Abstract 

We present a new method for estimating the frontier of a multidimensional sample. The 
estimator is based on a kernel regression on the power-transformed data. We assume that 
the exponent of the transformation goes to infinity while the bandwidth of the kernel goes 
to zero. We give conditions on these two parameters to obtain complete convergence and 
asymptotic normality. The good performance of the estimator is illustrated on some finite 
sample situations. 
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1 Introduction 

Let (Xi,Yi), i = 1, . . . , n be independent copies of a random pair (X, Y) with support S defined 
by 

S = {(x,y) £ExM;0<y<g(x)}. (1) 

The unknown function g : E — >• R is called the frontier. We address the problem of estimating g 
in the case E = R. d . Our estimator of the frontier is based on a kernel regression on the power- 
transformed data. More precisely, the estimator of g is defined for all x £ K d by 

/ n In \ 1 /P 

9n(x) - Up + 1)Y1 K h{* ~ Xi)YP / J2 K h (x - X t )\ , (2) 

where p = p n and h = h n are non random sequences such that h — > and p — > oo as n — > oo. 
This latter condition is the key so that the high power-transformed data "concentrate" along the 



frontier. We have also introduced Kh(t) — K(t/h)/h d where K is a probability density function 
(pdf) on M. d . In this context, h is called the window-width. 

From the practical point of view, note that, compared to the extreme value based estimators 
ITCT1 [TBI [T31 [T71 [T5] , projection estimators [32] or piecewise polynomial estimators [23 [UJ , this 
estimator does not require a partition of S and is thus not limited to bi-dimensional bounded 
supports. Moreover, it benefits from an explicit formulation which is not the case of estimators 
defined by optimization problems [T2J such as local polynomial estimators [30] HH1 121] ■ From 
the theoretical point of view, this estimator reveals to be completely convergent to g without 
assumption neither on the distribution of X nor on the distribution of Y given X — x (see 
Section [3]). Note however that (p + l) 1 ' p — > 1 when p —> oo. In fact, this correcting term is 
specially designed for the case where Y given X — x is uniformly distributed on [0,g(x)]. In this 
latter situation, the estimator is asymptotically Gaussian with the rate of convergence n~ Q /( d+a ) 
(see Section 3]). This rate is proved to be minimax optimal for a— Lipschitzian d— dimensional 
frontiers |26j . Chapter 5. This result is generalized in |27j to boundaries of more general regions. 
Other extensions are provided in jTSJ [UJ to densities of Y given X = x decreasing as a power of 
the distance from the boundary. We refer to [3 [3 [IT] for the estimation of frontier functions under 
monotonicity assumptions, and to [TJ [3] for the definition of robust estimators in this context. 
We conclude this paper by an illustration of the behavior of our estimator on some finite sample 
situations in Section [5] and by describing our future work in Section [6] Technical lemmas are 
postponed to the appendix. 

2 Notations and assumptions 

To motivate the estimator ([3]), consider the random variable Z = (p + 1)Y P and the conditional 
expectation r n (x) = K(Z\X = x). Estimating the frontier g is often related to estimating the 
regression function r n . For instance, if Y given X = x is uniformly distributed on [0, g(x)], we 
have r\j v (x) = g{x). A similar remark is done in |23j where regression estimators are modified to 
build estimators of the frontier, but the profound difference here is that p — > oo. This condition 
allows to obtain rl/ p (x) — > g(x) even when Y given X — x is not uniformly distributed (see 
Lemma [TJ below). We denote by / the pdf of the random vector X and we introduce 



where Zi = (p + Note that <fi n {x) can be seen as a classical kernel estimator of ip n (x) = 

f(x)r n (x) but keep in mind that p —¥ oo. Similarly, 




(3) 




(4) 



is an estimator of f(x) and 

r n (x) = 0n{x)/ f„(x) (5) 
is an estimator of r„(x). Collecting (J3j) , (0]) and J5]), our estimator ([2]) can be rewritten as 

g n {x)=r n {x) l 'v. 

To establish the asymptotic properties of g n {x), the following assumptions are considered: 
(A.l): g is a-Lipschitz, / is /3-Lipschitz, with < a < /3 < 1, 
(A.2): < g min < g(x), Vz e R d , 
(A.3): /(x) < / max < co, e M d , 

(A. 4): X is a Lipschitzian pdf on R d , with support included in B, the unit ball of R d . 
Note that (A. 4) implies that, for all q > 1, we have < J B K q {x)dx < +oo. 



3 Complete convergence 

In this section, the complete convergence of the frontier estimator toward the true frontier is 
established. The next lemma can be seen as the intuitive justification why no assumption on the 
conditional distribution of Y given X is required in the proof of Theorem [T] 

Lemma 1 Under (A.2), for all x G B, r n (x) 1 / p —> g(x) as n — » oo. 

Proof : Let e > 0. Since (X, Y) has support S defined by ([T]), it follows that 

r n {x) = (p + 1)1 (Y p \X =x)<(p+ l)g p (x) 

and thus, since (p + l) 1//p — s> 1 as p — >• oo, for n large enough and all x E B, 

r 1 J"(x)<(l+e)g(x). (6) 

Moreover, we have, 

r n {x) > (p+l)E(Y p l{Y>g{x)-e}\X = x) 

> (p+ l)(g(x) - e) p P(Y > g(x) - e\X = x). 

Now, since (X, Y) has support S, one can assume without loss of generality that Y given X = x 
has support [0, <7(x)] such that P(Y > g(x) — e\X = x) > 0. It follows that 

[(p + 1)P(Y > g{x) - e\X = x)] 1/p -> 1 

as p — > oo, and consequently, for n large enough, 

r 1 J p (x)>(l^e)g(x). (7) 

Collecting © and (J7J gives the result. ■ 



Theorem 1 Suppose (A.1)-(A.4) hold and nh d / \ogn — > oo. Then g n (x) converges completely 
to g(x) for all x £ M. d such that f(x) > 0. 

Proof : Let x eR d such that f(x) > and let e such that < e < g(x). Define < rj < 1/4 by 
7] = e / {Ag{x)). Then, from Lemma |51 

f r^ix) 



{\dn(x) -g(x)\ > e} 



- 1 



> 4r) 



C 



f(x)gP(x) 



> V } u 



fn(x) 



> V 



Since f n (x) converges completely to f(x), see e.g. [2], Chapter 4, Theorem III. 3, it follows that 
(/„(x)//(x)) 1 / p converges completely to 1. Therefore, writing 

i/p 



with 



T n {x) 



Y V 1 



i i/p 



g(x)J f(x) 



1 

-V K h (x-Xi 
n 

i—l 

and remarking that (p + — > 1 as n — > oo, it suffices to consider 

{|T n (a;) - 1| > n} C {T„(x) > 1 + r/} U {T„(x) < 1 - 77} . 
The two events are studied separately. First, let < 6 < rj. Then, \\x — < h entails 

Y t - g(x)(l +S)< g(Xi) - g(x) - Sg{x) < L g h a - Sg min < 
for n large enough and where L g is the Lipschitz constant associated to g. We thus have 



T n {x) 



1 " 

-Y^Khix-Xi) 



i=l 



!; ^l{Yi<g(x)(l + S)} ' 



< (1 + *) 



1 1 

- £ - < g(x)(l + 5)} — 



/(*) 
i/p 



and consequently, 



{T n (x)>l + V } C |l^^ a; -X i )l{F i < 5 (a ; )(l + 5)}-^> 

since, for n large enough, ((1 +ry)/(l + ^)) p > 2. From [2], Chapter 5, Corollary II. 4, the following 
complete convergence holds: 



1 - 1 
- J2 ~ XiMYi < g(x)(l + 5)} — 

o — 1 J \ ) 



\Y < g(x)(l + S)\X = x) = 1, 



and therefore 



oo oo / ^ n ^ 

HTn(x) >1+^)<E P -^ Kh{X ~ Xl)1[Yl K 9{X){1 + 6)} 7(X) 
n=l n=l \ i=l ■' ^ ' 



which concludes the first part of the proof. Second, 

v \P 

1 t 



T n (x) > 



1 " / Y \ p 

-^Knix-Xi)^) m>g( X )(l-S)} 



1 



1 > 1 < +oo, 



1 i/p 



\ ]T K h (x Xi)l{Yi > g(x)(l S)}-± 



i/p 



i=l 



and consequently, 



{T n (x) < 1 - ry} C 1 1 £ K h (x X t )l{Y > g(x)(l <J)}-L < (^j} ■ 
Now, since P(Y > g{x)(l — S)\X = x) > 0, there exists 7 > such that, for n large enough, 

' ~''^ P -F(Y>g(x)(l-5)\X = x)<-i, 



entailing that, for n large enough, 

{T n (x) < l-r]} C |i J2 R h( x - Xi)l{Yi > g(x)(l - 5)}-^ - ¥(Y/g(x) > 1 - S\X = x) < -7 
Taking into account of the following complete convergence 

n 1 

- ^ Jf fc (a; - Xi)l{Yi > g(x)(l - 5)}— -A P(y > - 6)\X = x), 

n i=1 /W 

it follows that 

00 

]TP(T„(x) < 1 - r?) 



n=l 



00 / n 

n=l \ i=l ^ ^ ' 



¥(Y > g(x)(l - 6)\X = x) < -7 



< +00, 



which concludes the second part of the proof. 



4 Asymptotic normality 

Second, the asymptotic normality of the frontier estimator centered on the true frontier is estab- 
lished. To this end, asymptotic expansions of the expectation and variance of ip n (x) are needed. 
These calculations are done under the additional assumption 
(A. 5): Y given X = x is uniformly distributed on [0, g(x)). 



The next two lemmas are similar to classical ones in kernel regression (see for instance [5], The- 
orem 6.11), but the dependence on n of the function ip n (x) induces technical difficulties. We 
first establish that <p n {x) is an asymptotically unbiased estimator of (p n (x) in the sense that 
Eip n (x) I ip n (x) — > 1 as n — > oc provided that ph a —> 0. 

Lemma 2 Under (A.1)-(A.5), if ph a -> 0, then for all x e R d 

E<p n (x) = <p n {x) [1 + 0(ph a )} . 

Proof : From ([3]), it follows that 

Etp n (x) = E(K h (x - X)Z) = E(K h (x - X)E(Z\X)), 
so that, by a straightforward calculation, and recalling that ip n (u) — g p (u)f{u), we obtain 

E(p n {x) = E(K h (x-X)gP(X))= J^^K (^^) <p n (u)du (8) 



B 



K(y)cp n (x - hy)dy, 

with a classical change of variable, and since K has a compact support. We thus can write: 

E<p n (x) - cp n (x) = / K(y)[<p n (x - hy) - <p n (x)]dy. 
Jb 

Consider now the decomposition below: 

\<p n (x - hy) - y n {x)\ < f{x - hy) \g^(x - hy) - g p (x)\ + g p (x) \f(x - hy) - f(x)\ := 7\ + T 2 . 
Following Lemma [SI 



T x = f(x~hy)g p (x) 



g p {x -hy) _^ 



< 2f Tnax ^-g p (x)ph a = g p (x)0(ph a ), 

9m'm 



gP(x) 

T 2 < 9 p {x)L f h? = g p (x)0(h^) = g p (x)o(ph a ), 
where Lf and L g are the Lipschitz constants of the functions / and g. Finally, 

E(p n [x) - <p n (x) = g p (x)0(ph a ) = <p n (x)0(ph a ), 

and the conclusion follows. ■ 

Similarly, we now provide an equivalent expression for Y((p n (x)/cp n (x)) which appears to be of 

order p/(nh d ). Thus, condition p/(nh d ) — > will be necessary in Corollary [T] to obtain the weak 

^ ^ p 

consistency of <p n (x), i.e. to ensure that <p n (x) / (p n (x) — > 1. 



Lemma 3 Under (A.1)-(A.5), if ph a -> then for all x e R d , 



Proof : We have 

1 " 1 
Y((p n (x)) = —Y / V(K h (x-X l )Z l ) = -Y(K h (x-X)Z) 

1 — 1 



' ( ^-A) Z 2 ) - E2 (£»(*)) := T 3 + T 4 . 



nh 2d \ \ h J J n 
From Lemma [H we immediately derive 

T A = ^ 2 n {x) [l + o(l)] 



Ti = W^ I K>( S ) ds 4£[l + o(l)], (9) 



We shall prove that 

To = , v , , , . , , 

nh d 2p + 1 J B w f(x r 

leading to T4/T3 = 0(h/p), and the announced result follows. To this end, remark that 
T3 = ^ E ( / ' 2 (^) E(Z2m 

1 fe+^/W^W) 



R2{s)ds L (nr) 9 2p ( u )f(^ 



nh 2d 2p+l \ \ h 
n/i d 2p+l J B 

where we have introduced the kernel Q — K 2 / j B K 2 (s)ds. It is easily seen that the second integral 
is similar to this appearing in K0 n (x), (see (|5|1), with K replaced by Q and p by 2p. Thus, as in 
the proof of Lemma [2] we have 



and © is proved. 

As a simple consequence of Lemma [5] and Lemma [3J we have 

Corollary 1 Under (A.1)-(A.5), if ph a ->• andp/(nh d ) -4 0, tfien, /or aZZ a; G K d , 

<p n (x)/(p n (x) -t 1. 

We can now turn to our main result. 

Theorem 2 Suppose that nph d+2a andp/{nh d ) — > 0. Le£ us define 
Then, under (A.1)-(A.5), for all x G R d , 

(If - 4 ^ 



Proof : First, note that nph d+2a — > and p/(nh d ) — > imply ph a — > 0. From Lemma GO it 
suffices to prove that 



To this end, define 



(T~ 1 (x) 1 

Wi, n = " K ' , , K h (x - Xt)Zi 
np tp n (x) 



so that we can write 

n 
i=l 

Following Lemma [31 we have 



V(£ n ) = nY(W hn )= " 2 V ; 2 ^ V(y n (x)) 



(2p+l)nh d f(x) 1 f (p + lf f 2 y n (x) 

p*S B K*(s)ds vl(x)nh d 2p + f J g [S)dS f(x) + 
= l + o(l). 

Thus, the condition of Lyapounov reduces to 

n 

E \Wi, n - EW^ n \ 3 = nE \W x , n - EW hn \ 3 -»• 0. (10) 

Taking into account that W\ tTl is a positive random variable, the triangular inequality together 
with Jensen's inequality yield 

EIW^-EW^ 3 <8E(Wi 3 in ). 
Introducing the kernel K 3 / J„ K 3 (s)ds, and mimicking the proof of Lemma[3l we obtain 

E(WU = h »* (1 + 0(1)) = Kn-W^l + o ( l)), (11) 

where k is a positive constant. Returning to (ITU1) . we have 

^ £ |W 4 ,„ - EWi, n | 3 < 8k 1/2 (1 + o(l)) -> 0, 

i=l 

and the result is proved. ■ 
Remark 1 Theorem^ holds when a~ 1 (x) is replaced with 

since m this context f n (x) — > f(x). This allows to produce pointwise confident intervals for the 
frontier. 



Remark 2 To fulfill the assumptions of Theorem^ one can choose h = n 1 /( d + Q ) a nd p 
e n n a ^ d+a \ where (e„) is a sequence tending to zero arbitrarily slowly. These choices yield 




B 




which is the optimal speed (up to the e n factor) for estimating a— Lipschitzian d— dimensional 
frontiers, see f26f . Chapter 5. 

The good performances of g n (x) on finite sample situations are illustrated in the next section. 
Remark [2] will be of great help to choose p and h sequences. 

5 Numerical experiments 

Here, we limit ourselves to unidimensional random variables X (p = 1) with compact support 
E = [0, 1]. Besides, Y given X = x is distributed on [0, g(x)] such that 



with 7 > 0. This conditional survival distribution function belongs to the Weibull domain of 
attraction, with extreme value index —7, see [6] for a review on this topic. The case 7 = 1 
corresponds to the situation where Y given X = x is uniformly distributed on [0, g(x)]. The larger 
7 is, the smaller the probability (TT21 is, when y is close to the frontier g(x). The behavior of the 
proposed frontier estimator is investigated on different situations: 

• Two distributions are considered for X: a uniform distribution U([0, 1]) and a beta distribu- 
tion S(2,2). 

• Two frontiers are introduced. The first one 




(12) 



01 0*0 



1 + exp (-60(x - 1/4) 2 ) if < x < 1 /3, 

1 + exp (-5/12) if 1/3 < x < 2/3, 

1 + 5 exp (-5/12) - 6 exp (-5/12)x if 2/3 < x < 5/6, 

6x-4 if5/6<x<l. 



is continuous but is not derivable at x = 1/3, x 



2/3 and x 



5/6. The second one 



g 2 ( x ) = (1/10 + sin(7rir)) (11/10 - exp (-64(x - 1/2) 2 ) /2) 



is C 



• Four sample sizes are simulated n g {200, 300, 500, 1000}. 



• Three exponents are used 7 £ {1,2,3}. 



The following kernel is chosen 

K{t) = cos 2 (7rf/2)l{i g [-1, 1]}, 

with associated window width h = 4:a{X)n~ 1 / 2 and with p — n 1 / 2 . The dependence of these 
sequences with respect to n is chosen according to Remark [5] with a = d = 1 . The multiplicative 
constant Aa{X) in h is chosen heuristically. The dependence with respect to the standard-deviation 
of X is inspired from the density estimation case. The scale factor 4 was chosen on the basis of 
intensive simulations. 

Here, the experiment involves several steps: 

• First, m — 500 replications of the sample are simulated. 

• For each of the m previous set of points, the frontier estimator g n is computed. 

• The m associated L\ distances to g are evaluated on a grid. 

• The mean, smallest and largest L\ errors are recorded. 

Some results are depicted in |16j . where the best situation (i.e. the estimation corresponding to 
the smallest L\ error) and the worst situation (i.e. the estimation corresponding to the largest L\ 
error) are represented. Note that, even in the worst situations, the empirical choices of sequences 
h and p seem satisfying for all the considered frontiers and densities of X. In fact, the worst 
situations are obtained when no points were simulated at the boundaries of the support. 
Finally, the above estimator is compared to three other ones: 

• The estimator g n with p = 1, which reduces to a rescaling of the regression estimator, in a 
similar spirit as in [23 . 

• Geffroy's estimator [10], denoted by g^, which is a step function based on the extreme values 
of the sample. 

• The kernel estimator g^ introduced in |15j . which is a smoothed and bias-corrected version 
of Geffroy's estimator. 

Results are summarized in Table [TJ It appears that, when 7 increases, performances of all esti- 
mators decrease, since the simulated points arc getting more and more distant from the frontier 
function. In the case p = 1 and 7 = 3, one can see that g n does not converge to the true frontier 
when n increases. This shows that the condition p — > 00 is necessary to obtain the convergence 
of the estimator. Finally, note that in all the situations considered in Table [TJ g n outperforms g^ 
and g%. 



6 Conclusion and further work 



To conclude, let us note that, even though g n converges to the true frontier g in case of non uniform 
conditional distributions, it is possible to design new estimators dedicated to particular parametric 
models. For instance, in case of model (fT2|) . estimator g n could be modified to obtain 

^ W= ( 7 B(1+M) t,KH(x-Xt)Y? I , 
where B is the beta function defined by 

B(a,b) = f u^il-uf^du. 
Jo 

Of course, g n corresponds to the particular case g n ,i- When 7 is assumed to be known, the 
new multiplicative constant yields a very efficient bias correction, see |16j for an illustration. A 
part of our future work will consist in defining an estimator of 7 and plugging it into g n ,7- New 
asymptotic results will be established. We also plan to investigate the asymptotic properties of 
local polynomial estimators based on the same ideas as those used for g n and <?„ i7 . 



7 Appendix: Auxiliary lemmas 

The following lemma provides convenient bounds obtained by a specific study of the functions 
u — > |(1 + u) p — 1| — 2p\u\ and u — y (1 + u) 1 / p — 1 — ^u. The study is left to the reader. Note that 
these bounds could not be directly derived from the Taylor formulas |(1 + u) p — 1| = \pu + o(u)\ 
and (1 + it) 1 /? — 1 — In — ~ I)" 2 + °( u2 ) where the dependence on p of o(u) and o(u 2 ) is 

not precised. 

Lemma 4 Suppose p > 1. 

(i) Then, p\u\ < In 2 entails (1 + u) p - 1 < 2p\u\. 



(ii) Let C>2. Then, \u\ < 1/2 entails (1 + u) 1/p - 1 - -u 



The next lemma is dedicated to the control of the local variations of the frontier on a neighborhood 
of size h. 

Lemma 5 Suppose (A.l), (A. 2) hold. If ph a —> and \\x — y\\ < h, then for sufficiently large n, 



3(-)^ p _ 1 



My) 

where L g is the Lipschitz constant of the function g 



9 'min 



Proof : Take u — — 1 and observe that p\u\ < Pj^-\\x — y\\ a , Thus, if \\x — y\\ < h, and 
ph a —¥ 0, we have p\u\ < ln2 for sufficiently large n. Then, following Lemma 0Ji) , for sufficiently 
large n, we obtain 



(l + u)P-l\ 



ffO) 

g(y) 



< 2p\u\ < 2—?-ph a , 



and the result is proved. ■ 
Lemma [6] is used to establish the complete convergence of random variables ratio. 

Lemma 6 Let S, T be real random variables, a, b non zero real numbers, and < rj < 1/2. Then, 



S a 
T ~ 1 



> Arj 



C 



S --l 

a 



T 

J' 1 



Proof : Consider the following obvious equality: 



= 11-? 



S a 
T~b 



> V 



T 



(13) 



The triangular inequality yields for all 77 > 0: 



S 

1 

a 



<r]}n 



T 

J' 1 



S a 
T ~ 1 



< 2n 



Taking < 77 < 1, we obtain 



a 



< v ? n 



T 

J' 1 



S a 
T~b 



< 



2r? 
1-rj 



S a 
T ~ b 



Finally, note that < Ar\ for < n < 1/2. 



The next three lemmas are of great use to deduce successively the asymptotic normality of g n (x) 
from r n (x) and the asymptotic normality of r n (x) from !p n (x). 



Lemma 7 Let x € 



If fn(x)/ f(x) 1 and <p n (x) / '<p n (x) -> 1, then 



r n (x) 
r n {x) 



-1 = 



&n(x) 
¥>n(x) 



- 1 



/(*) 



-1 (l+Op(l)). 



Proof : The hypotheses yield = f^/7^7 4 



1. Thus it suffices to consider S = Vn< ) x \ 



T = 



/GO 



, and a = b = 1 in the equality (|13l) . 



Lemma 8 Let x € R d . If f n (x)/f(x) 1 erne? (p n (x) / (p n (x) — > 1, £/ien 



-1 =- 



1 / r n (») 



P V 7 ""^) 



-1 (l + o p (l)). 



Proof : From the hypotheses, w n (x) := ffifj — 1 = o p (l). Moreover, following LemmalUJii), on 
the event {|to n (x)| < 1/2} we have: 



A n {x) 



9n(x) 
9{x) 



- 1 )- - ( ^4 - A - a + -«(-)) 1/p - 1 - — 

/ P V^nW / P 



P 



We thus have, on the one hand, 



pA n (a;)l{| Wn ( a: )| < i/2} = o p (w n (x)). 
On the other hand, for all e > 0, 

L^M 1{lwn(x) ^ 1/2} > £ | c {\w n (x)\ > 1/2} 

leading to 

P (p^7^1{| W „(x)|>l/2} > 4 < P{\Wn(x)\ > 1/2} -> 0, 

and thus 

P 7-yl{|t«„(x)|>i/2} = o p {w n {x)), 

which completes the proof. 

Lemma 9 Suppose that nph d+2a — > and p/(nh d ) — > 0. Let us define 

1/2 



*'M-«* +1 >^(7^j*) • 



and let Q be an arbitrary distribution. Then, under (A.1)-(A.5), 



p \Vn{x) <Pn(x) J J I " V 9( x ) 

Proof : First, note that nph d+2a —> and p/(nh d ) — > imply pft" — > 0. Thus, from Corollary [TJ 
tp n (x) / (p n (x) -H- 1. Besides, p/(nh d ) — > implies n/i d — > oc, and thus, using a classical result 
on density estimation (see for instance [2j, Chapter 4, Theorem II. 1), we have f n (x)/f(x) —> 1. 
Lemma [7] thus entails 

gV*lg) ^n(x) \ CTrT 1 ^) Z' gn(g) ,\ f /~M M V> 

^Um _1 J = — U^rV — —{-m- 1 ) il + 0p{1)) 

gn^) f _ E^n(fO \ _ C^Q) / /n(x) _ Ef n {x) 

p \tp n (x) (p n (x) J P \ /(- T ) 



(1 + 0,(1)) 



Following Lemma [H we have, 

^(x) ^(») ^ =0 ff!^y /2> \ (p 0=0 (( np ^) 1 / 2 ) =o(i), 



P V </>n0) / \\ V ) J 

and from a classical result on density estimation E/ n (x) — /(x) = 0(h a ), see [4], Proposition 2.1, 
we have 

a-Hx) f^)_ 1 ) =0 ((^) 1/2 ) 0(ha) = r {n ^^ 



p \ / v V p 



Consequently, 



P \r n (x) J p 



(x)f? n ( X ) 1 \_^MIfn^)__ 1 \ {1 + 0p{1))+0p{1) 



Vn(x) J p \ f(x) 



Again, using a classical result on density estimation, V(/„(x)) = 0(l/(nh d )), see g], Proposi- 
tion 2.2, we have 



a-'(x) (r„(x) _\ = o^) / vJ^l _ j | + (14) 



and thus 

P \r n (x) ' ) p \ip n (x) 

Suppose now that there exists a probability distribution Q such that 

Vn 1 (x) f <Pn{x) E(fi n (x)\ ^ 



p \<Pn{x) <Pn(x) 

From (fl"4)) . we deduce that 

VnHx) f r n {x) ^ 4 Q 

p V^W / 

Finally, from Lemma [8] we can conclude that 



! , , I 9n(x) A j „ 



and the result is proved. 
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7=1 


n 


g n with p — > oo 


g n with p = 1 


9* 


9n 


200 


0.121 [0.051, 0.237] 


0.651 [0.407, 0.907] 


0.134 /0.050, 0.261] 


0.183 /0.050, 0.55^/ 


300 


0.100 /O.O.^, oj5^/ 


0.636 /0.^5, 0.831] 


0.111 /O.OOi, 0.^9/ 


0.157 /0.075 ; 0.300] 


500 


0.078 [0.042, 0.138] 


0.627 /O.^i, 0.5i5/ 


0.087 [0.046, 0.168] 


0.128 /0.00^ 7 0.234] 


1000 


0.057 /0.025, 0.112] 


0.616 [0.486, 0.752] 


0.062 /0.055, 0Ji7/ 


0.093 /O.O^O, 0J55/ 


7 = 2 


n 


g n with p — > 00 


g n with p = 1 


9« 




200 


0.321 /0.107, 0.^30/ 


0.575 [0.415, 0.759] 


0.337 /0J50, 0.5i0/ 


0.426 [0.269, 0.591] 


300 


0.297 [0.194, 0.457] 


0.562 /0.500 ; 0.755/ 


0.311 [0.171, 0.490] 


0.393 /0.255, 0.569] 


500 


0.262 /OJtffil, 0.379] 


0.545 [0.429, 0.667] 


0.275 /0J72, 0.550/ 


0.347 /0.25i ; 0.^52/ 


1000 


0.226 /0J55, 0.505/ 


0.533 [0.463, 0.623] 


0.240 [0.152, 0.336] 


0.293 [0.200, 0.388] 


7 = 3 


n 


g„ with p — > 00 


g n with p = 1 


9* 


9n 


200 


0.526 [0.331, 0.709] 


0.749 [0.627, 0.888] 


0.550 [0.340, 0.724] 


0.624 /O.^iO, 0.750/ 


300 


0.496 |0.565 0.669] 


0.744 /0.052, 0.865] 


0.523 /0.57i, 0.057/ 


0.591 [0.452, 0.739] 


500 


0.457 |0.506, 0.550/ 


0.741 /O.O^O, 0.5i7/ 


0.486 /0.575, 0.620] 


0.545 [0.434, 0.668] 


1000 


0.410 [0.315, 0.505] 


0.742 [0.685, 0.817] 


0.442 /0.527, 0.531] 


0.486 /0.575 ; 0.575/ 



Table 1: Comparison between L\ errors obtained with four different estimators. The mean error 
is indicated as well as the range between the minimum and the maximum error. The experiments 
are conducted on a 5(2, 2) covariate, with frontier function g 2 . 



